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Abstract 

Hamiltonian models based on two different infrared approximations are studied in order 
to obtain an explicit comparison with the standard analysis of the infrared contributions, 
occurring in the relativistically covariant perturbative formulation of Quantum Electrody- 
namics. 

Moller operators, preserving respectively the Hilbert scalar product, for the Coulomb-gauge 
models, and an indefinite metric, for the models formulated in Feynman's gauge, are ob- 
tained in the presence of an infrared cutoff, after the removal of an adiabatic switching and 
with the aid of a suitable mass renormalization. 

In the presence of a dipole approximation, spurious contributions to the infrared factors 
are shown to necessarily arise in Feynman's gauge, with respect both to the Coulomb-gauge 
model and to the amplitudes of Quantum Electrodynamics, and the connection of this result 
with a recent work on the Gupta-Bleuler formulation of non-relativistic models is discussed. 
It is finally proven that by dropping the dipole approximation and adopting an expansion 
around a fixed charged particle four-momentum, first introduced and employed in the study 
of the infrared problem by Bloch and Nordsieck, the infrared diagrammatic is fully repro- 
duced and spurious low-energy effects are avoided. 

Key words: quantum electrodynamics; infrared problem; local and covariant gauge; indefi- 
nite metric; solvable models 
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Introduction 



In Quantum Electrodynamics {QED), the description of states at asymptotic times and 
the derivation of the scattering matrix are still open issues; it is customary to refer to such 
questions as the infrared ( IR ) problem. 

At the perturbative level, transition amplitudes between states containing a finite number 
of photons are ill-defined, since radiative corrections due to soft photons typically exhibit 
logarithmic divergences Wci95, JR76 . As a consequence, in contrast with ordinary quantum 
field theories, Dyson's S - matrix |Dys49| |Dys51| is defined only in the presence of a low- 
energy cutoff and the problem of a proper identification of asymptotic states arises. 

As early as 1937, in their pioneering paper on the subject BN37 , Bloch and Nordsieck 
proved that IR singularities arise in the perturbative expansion because of basic physical 
facts; they argued that, on the basis of the correspondence principle, one has to expect a 
vanishing probability for the emission of a finite number of photons in any collision process 
involving electrically charged particles. 

Exponentiation of the low-energy photon contributions was conjectured by Schwinger 
[Sch49j and then proved by Yennie, Frautschi and Suura (YFS) in [YFS61] within the 
framework of the local and covariant Gupta-Bleuler formulation |Gup50| QBle50 of QED . 
This led to a pragmatic approach to circumvent the soft-photon divergences; one introduces 
an IR cutoff, sums the transition rates over all final photon states with energy below the 
threshold of the photon detectors and removes the regularization in the resulting expression. 
The finiteness of the so-obtained inclusive cross-sections is ensured by cancellations, at each 
order in perturbation theory, between the virtual IR divergences and those due to soft- 
photon emission Wei95l IJR76] . 

It is important to remark that such a recipe somehow avoids to take into account the 
properties arising in the characterization of physical charged states, such as the spontaneous 
breaking of the Lorentz boosts in the charged superselection sectors and the absence of a 
sharp eigenvalue for the mass operator of an electrically charged particle. These features 
were established in past decades through several mo del- independent investigations [SW74, 
IFPS741 IFPS771 IFMS79al IFMS79bl IBuc82l IBuc86j . mostly within the algebraic approach to 
quantum field theory and independently of the perturbative-theoretic framework. 

Quite generally, since its early days perturbation theory has been mostly applied to 
extract physical predictions (notable exceptions are |Sym71[ IMS 8 3 , Stcinb] ) and its relation 
with the above-mentioned structural (non-perturbative) properties of QED is still unclear; 
in particular, local and covariant quantizations are incompatible with positivity |St67| and 
are only consistent with a generalization of the Wightman axioms. 

The above discussion provides motivations to understand and possibly fill the gap that 
separates the Feynman-Dyson formulation from a collision theory in which the structural 
aspects of the infrared problem are taken into account; such an issue is important both 
conceptually and practically, since perturbation theory remains the only source of detailed 
information on QED and its local and covariant version is the best controlled one regarding 
the renormalization procedure. 

As a starting point to address this issue, in the present paper we seek for a mathemati- 
cal formalization of the hypotheses underlying the local and covariant treatment of the IR 
divergences occurring in QED . In particular, we compare the infrared amplitudes of the 
standard analysis with those obtained from solvable hamiltonian models based on two dif- 
ferent approximations, which might seem to be equally suited, from a physical point of view, 
for an investigation of low-energy photon effects; the electric dipole approximation and the 
expansion around a fixed (asymptotic) charged particle four-momentum. 

Non-relativistic Coulomb-gauge models with a dipole approximation have been studied 
rigorously in Hl;i(i9. Ara83 in order to obtain non-perturbative constructions, in particu- 
lar of asymptotic electromagnetic fields and of one-particle charged states. More recently 
[HS09 , asymptotic e.m. fields and a scattering operator have been constructed within the 
Gupta-Bleuler formulation of non-relativistic QED , still in the presence of a dipole approx- 
imation, by addressing the questions connected with the existence and uniqueness of the 
dynamics of Heisenberg operators in the presence of an indefinite metric. 
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The problem of the comparison of such models with the perturbative structures and 
results has however not been addressed explicitly, thus leaving the validity of perturba- 
tive procedures in a dubious condition. In addition, as we shall see, the Gupta-Bleuler 
formulation of non-relativistic models with a dipole approximation suffers from substantial 
limitations in the comparison with the Feynman-Dyson expansion. 

Our main result |PhDth] is that (only) within hamiltonian models of the "four-vector 
Bloch-Nordsieck type" the infrared diagrammatic of the local and covariant expansion of 
QED is reproduced by Mdller operators, obtained as weak limits in the presence of a low- 
energy regularization and a suitable mass counterterm, after the removal of an adiabatic 
switching. In contrast, the dipole approximation prevents to recover fundamental features 
of the infrared contributions to Feynman's amplitudes; basically, its effects on current con- 
servation, also recognized and studied in [HS09 , give rise to substantial additional contri- 
butions to the Moller operators, resulting in particular in a discrepancy by a factor 3/2 
with respect to the Coulomb-gauge result (and to covariant perturbation theory in QED ) , 
in the exponential factors describing soft-photon effects. 

The structure of the paper is as follows. Section[TJis devoted to the analysis of the Pauli- 
Fierz model |FP38) . a non-relativistic model formulated in the Coulomb gauge and based 
on the dipole approximation, taking Blanchard's treatment [Bla69l as a starting point. We 
prove the existence and uniqueness of the dynamics and obtain Moller operators as strong 
limits of the corresponding evolution operator in the interaction representation, for a fixed 
value of an infrared cutoff. 

In order to set up a comparison with the perturbative procedures and results we in- 
troduce a four- vector model, retaining the approximations of the Pauli-Fierz Hamiltonian. 
Within a simpler setting with respect to [HS09 , we define the dynamics of the model and 
prove its uniqueness; then we construct Moller operators as weak asymptotic limits of the 
evolution operator (in the interaction representation) for a fixed infrared cutoff, and discuss 
the spurious effects induced by the dipole approximation. 

In Section[5]we introduce models based on an expansion already implicit in |BN37j . here- 
after referred to as Bloch-Nordsieck (BN) models. We will first consider a model formulated 
in the Coulomb gauge and then its four-vector version. The existence and uniqueness of the 
dynamics of the models and the control of the asymptotic limits of the evolution operators 
are achieved within the same framework of Section [TJ We show that the infrared amplitudes 
of QED are reproduced by means of Moller operators and that the contributions arising in 
the presence of the dipole approximation disappear. 

We conclude the paper with an outlook for future research. 

Notations 

The metric g p v — diag (1, — 1, — 1, — 1 ) of Minkowski space is adopted and natural units 
are used (fi = c=l).A four- vector is indicated with v ^ or simply with v , while the 
symbol v denotes a three-vector. We use the symbol c • d for the indefinite inner product 
between four- vectors c and d . 

The norm of a vector E L 2 is indicated by || (j> || 2 , the Hilbert scalar product by ( . , . ) 
and indefinite inner products by (.,.). 

The Hilbcrt-spacc adjoint of an operator .4 is denoted by A *, while the symbol B t stands 
for the hermitian conjugate, with respect to the indefinite inner product (.,.), of an oper- 
ator B defined on an indefinite-metric space. 

We denote by & the symmetric Fock space, by | . | the norm of & and by N the number 
operator. The projection of <f> € & onto the n - particle space is denoted by <fi ^ = S n <f> , 
where the orthogonal projection S n is the symmetrization operator defined in terms of the 
permutation group of degree n . 

In the Coulomb-gauge formulation, a s ( a * ) will stand for the photon annihilation (cre- 
ation) operator-valued distribution, fulfilling the canonical commutation relations (CCR) 

[o,(k), o:,(k')] = 6 SS , «(k-k'), a) 
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with s and s ' polarization indices. 

In the same gauge, the Hamiltonian of the free e.m. field is denoted by i?o'™' and the 
vector potential at time t — by 

A tr (x)= Yl I -TW^= e ^ k ) t a ^ k ) e ' ik ' X + <( k ) e ~ lk ' X ], ( 2 ) 

where e s , s = 1, 2, are orthonormal vectors, satisfying the transversality condition k • 
e s (k)= 0. _ 

The annihilation and creation operator-valued distributions in the FGB gauge, denoted 
respectively by a» ( k ) and a ^ ( k ) , fulfill the CCR 

[a"(k), a"t( k ')] = -gK" S(k-k'). (3) 

In the same gauge, the Hamiltonian of the free e.m. field is denoted by H § m - and the 
vector potential at time t = by 

A"(x) = / d3fc [o"(k) e jkx + o"t(k) e- 1 kx ] . (4) 

J V 2 W k 

We write 

F""(x) = 9^ A ,y (x)- 5 l/ ^l^(x) (5) 

for the e.m field tensor at t = . 

The convolution with a form factor p is indicated by 

x) = | d 3 C P (0 A"(x-f), (6) 

and similarly for A tr . For brevity we write 

a(f(t)) ee J d 3 k /*(k) / p (k, (7) 

and denote the corresponding sum in the Coulomb gauge by a tr (/(<)) ■ 

( R 3 ) will stand for the Schwartz space of C 00 functions of rapid decrease on R 3 . 



1 Pauli-Fierz-Blanchard Models 

In the present Section we discuss the Pauli-Fierz-Blanchard (PFB) model and formulate a 
suitable four-vector model, retaining the approximation of the Pauli-Fierz hamiltonian. 

The model |FP38j describes the interaction of a spinless Schrodinger particle with the 
quantum electromagnetic field, under suitable infrared approximations. It was reconsidered 
three decades later by Blanchard, who investigated the questions connected with a mathe- 
matical formulation of the fact that an infinite number of photons is emitted in any collision 
process involving electrically charged particles. In [Bla69l , he proved the existence of the 
dynamics and showed that a unitary operator can be obtained as the limit of evolution op- 
erators in the sense of morphisms of a suitable (C*-) algebra. Furthermore, he established 
the existence of Moller operators interpolating between the Pauli-Fierz Hamiltonian and its 
perturbation by a potential term, for a large class of potentials. 

The model was further studied by Arai |Ara83) . who was able to construct one charged- 
particle states and obtained a scattering operator at fixed charged particle momentum, also 
allowing for the inclusion in the hamiltonian of the bilinear term in the gauge field. In the 
sequel we shall take Blanchard's treatment as a starting point, although our analysis will 
require some changes with respect to his setting. 

First, since we wish to employ the model in order to investigate the methods at the basis 
of the local and covariant perturbative treatment of the infrared divergences, we adopt an 
IR cutoff throughout the analysis. Even in its presence, the limits considered by Blanchard 
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for the Moller operators only exist in the weak topology and do not define unitary operators; 
in order to obtain strong convergence and unitarity, a mass renormalization and an adiabatic 
switching of the interaction will be essential. 

Second, we introduce a four-vector version of the model, prove the existence and unique- 
ness of the dynamics and control the large-time limits of the time-evolution operator in the 
interaction representation, yielding Moller operators as isometries on an indefinite vector 
space. 

We consider the infrared-regularized PFB Hamiltonian 

(PFB) P 2 

H in t,tr=-— P ' A tr(P> X=0) . (9) 
TO 

The particle, of mass to , charge e and rotationally invariant distribution of charge p G 
J? ( K 3 ) , will also be called electron. The subscript A on the left-hand side of © denotes 
the fictitious photon mass employed as an infrared-regularization method as in QED , by 
setting w£ = k 2 + A 2 , and the A - dependence in H g ™ r and in ([9]) is understood. The 
functional form of the interaction is dictated by the electric dipole approximation and implies 
that the electron momentum is conserved, while the total one is not. 

The Hilbert space of states of the model is ffl = L 2 ( R 3 )<g)jF , with L 2 the one-particle 
space and & the Fock space of photons. Since H is the sum of two positive and commuting 
self-adjoint operators, it is self-adjoint; in particular, it is essentially self-adjoint (e.s.a.) on 
D =\y{R 3 )(g>D Fo , with D F = (tp G F ; V (n) € S n <g> £ = 1 ^(M 3 ), Vn), F the 
set of finite particle vectors of & . 

Next we prove the essential self-adjointness of the Hamiltonian ((5J), making use of tech- 
niques already exploited in |Nel641 IAra81j . 

By the estimate |j b s (f) * || < || / 1| 2 \\(N + I) 1/2 * || , V f G Dom{N 1 / 2 ), where 
b s ( / ) j / £ L 2 , stands for either a s ( / ) or a * ( / ) , one finds 

||AUp,x = 0) *|| < c(p) WiN+l) 1 ' 2 *|| , VfeD fo , (10) 

with c(p) is a (positive) constant, provided the form factor p is held fixed. Further, for 
any A > 0, || (N+ I) *|| < A _1 || H e ;'™- * || + || * || , V* e D F , hence, by virtue of the 
estimate 

2||(A(g)B)$||<||(A 2 «)I + I«)B 2 )$||,V$eL» , 
with A = p,B=(A^+I) 1 / 2 ,one gets the bound 

\\Hi PFB) $|| < d(e, A;p) ||(ir +I)*||, V$€D , (11) 

for a suitable function d ( e , A ; p ) . By the bound (fTTj) and the CCR, there exists g ( e , A ; p) 
such that |($, [ffj PFB) ,iJ ] $)| < g(e, A; p)||(_ff + I) 1 / 2 $ || , V$ € J ; Nelson's 
commutator theorem, in the formulation given by Faris and Lavine |FL74[[RSII) . then implies 
that H ^ PFB '> i s e.s.a. on D , for all | e | , A > . 

Since the Hamiltonian commutes with the electron momentum operator, it is e.s.aj on 
(almost) any of the subspaces, on which p takes a constant value, obtained by decomposing 
D q on the joint spectrum of the components of p ; without losing generality we consider 
therefore p fixed in the rest of the analysis of the PFB model. 

The time-evolution in the interaction picture, corresponding to an Hamiltonian H = 
H + H int , with H o the free part and H int a time-independent interaction, is governed 
by the operator 

U j (t) = cxp(iH t) exp{-iHt) (12) 

and satisfies, on D , 

i = Hj(t) U^t), (13) 

The e.s.a. of the hamiltonian at fixed p also follows by the Kato-Rellich theorem. 
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with 

Hj(t)= exp (iii t ) H mt exp(-iH t). 

(|13p identifies U j (t) (a more general result, also applying to time-dependent interactions, 
is stated below). Since the commutator of H 1 evaluated at different times is a multiple of 
the identity operator at each definite momentum, we make use of the formula [Wil67 

e A + B = e A e B e -i [A,B] (w) 

Here and in the following, (fT4|) will be used to compute evolution operators, with an a 
posteriori verification that the so-obtained operators satisfy the evolution equation on a 
suitable domain, here on a core of the common domain of the operators H j (t) . We obtain 

U j (t) = exp (- i f dt' H z (*')) exp {-\ [ dt' 

■ t ' 



o 2 j o 



dt" [Hj (t'),Hj (*")]), (15) 

hence U (t) = exp ( — i H t) U j (t) . By explicit calculations one gets 

Uj(t)= c(t) eX p(ie(a* tr (f p (t)) +a tr (J p (t)))), (16) 

, i e 2 P 2 f d 3 k „ 9 . , N , sin u; t i N N , . 

c (t) = exp ( — -f- / — p 2 (k) (t —*-)), (17) 



3 to 2 J w k W k 

It is easy to check that \] I (t) does not converge for large times; thus we introduce a 
mass counterterm and a regularization of the oscillating terms occurring in (|17[) . (|18[) . by 
replacing the electric coupling by e ( ad -* ( t ) = e e ~ e ' * ' . The resulting Hamiltonian is 

tt(PFB) _ P 2 rj-e.TO. , „(e) a. Q s, 

H \,R - °< tr + H int,tr,R ' (- i9 ) 

„ 2 

n int , tr , R ~ n int , tr e ' z e 2 m ' 

4^ P 2 (k). (21) 



3 771 j lu k 

In the sequel, we state the results for positive times; in order to obtain the corresponding 
expressions for t < , it suffices to replace e by — e . Although the above Hamiltonian 
is time dependent, existence and uniqueness of the time-evolution unitary operators follow 
from the independence of time of their selfadjointness domain and strong differentiability of 
(if ( t ) - i)(H(0)- i)' 1 , through the results of [Kat56j . 
Upon inserting f|20[) . (I2l1) into the right-hand side (r.h.s.) of (|T5|) we obtain 

U$'l >x (t)= c< e >(t) exp (ze (a* tr (f^ (t)) + a tr (T p ie) (t)))) , (22) 

with 

,2p2 „2 „-2et 

3m 2 \ j j f \ 2 m 2 e 



c«"(() s exp(ii-E- <!«>(*)) exp(ie 2 j -E- -), (23) 



""'"'-/ Sf^ 1 '"'^''^ 1 '""- 11 '' <24> 



x/ 7 



k 
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(|23p .([M ]) provide a regularization of (fTTj) ; in particular, the oscillating term on the r.h.s. 
of (fM|) vanishes for t — > oo , due to the presence of the adiabatic factor, and the residual 
contribution of order 1 / e from d ( e ) ( t ) to the first exponential on the r.h.s of (|2U)) is can- 
celed, for e — > , by the z - dependent exponential. Hence the existence of the asymptotic 
time limits and of the adiabatic limit of (|23[) is proven. 

The strong convergence of the evolution operator (f2"2"j) can now be established: 

Proposition 1. By choosing the coefficient of the mass counterterm as in \21\) . both the 
large-time limits and the adiabatic limit of the evolution operator 122(1 . defining the Moller 
operators, exist in the strong topology of Jfi? : 

n±,*r = * ~ I™ lim U\% <x {-t) = 

e — » t — > =p oo ' ' 

= exp(- ie K (/ P .) + o.(7 P .)]). (26) 

5 



/ ps (k) = L : 



= L 



P(k) 



(0 

P ' £ s (k) 



- lim lim /^(k, t) 



//), 



— • (27) 



Proof. The limits of /ps (k, t) exist pointwise and thus in the strong L 2 topology by 
the dominated convergence theorem. 

The operator $(/ p e) (k, *)) = a (r (/ p e) ( k , t)) + a* r (/ p e) (k, t)) admits F as a 
dense and invariant set of analytic vectors and is therefore e.s.a. on F due to Nelson's ana- 
lytic vector theorem |RSIIj . Linearity and Fock-space estimates imply that $ (/p 6 " 1 (k, t)) 
converges strongly to 3>(/ p e '(k)) on F . Since <I>(/p e '(k)) is e.s.a. on F , one has 
convergence in the strong generalized sense and by Trotter's theorem |RSIj the existence of 
the time-limits in (f2l))) follows. A similar proof allows to establish the strong convergence of 
$(/^ e) (k)) for e^O. □ 

We wish to remark that without the adiabatic regularization the time limits of the operators 
considered above would only exist as weak limits; as a matter of fact, the Riemann-Lebesgue 
lemma implies the existence of the weak limits w — lim t^-^oo /ps~°^(k, t) = / ( k ) , 

while || /p S -0 ' (t) — f || 2 does not converge to zero. 

Treating the model as a description of the soft-photon effects in QED and introducing a 
unitary operator W on L 2 (R 3 ) , interpreted as a scattering operator for electrons, acting 
at time zero, a non-trivial S - matrix is given, at fixed A , by 

= s - lim lim P};| r , x (t)f < A (t'). (28) 

€ — > t , t ' — > +00 ' ' 

In order to compare the infrared diagrammatic of QED with the expansion of Moller's 
operators in powers of the electric charge, it is necessary to formulate a model in a gauge 
employing four independent photon degrees of freedom, such as Feynman's gauge. With 
this aim, we introduce, as in [H.S09 , a four- vector model retaining the approximations of 
the Pauli-Fierz Hamiltonian. 

The model is defined in the tensor product space "V = L 2 ( R 3 ) ® S? , Sfan indefinite- 
metric photon space to be constructed below. The Hamiltonian is given, on the domain Y 
defined in (1551) . by 

H {PFBR) = m + 1 m v 2 + jje.rn. _|_ g £ . A (p , X = 0) = H + H mt , (29) 

with w A ' = (l+v 2 /2,v),v = p/m, and will be referred to as PFBR Hamiltonian. 

In general, for models employing covariant vector potentials, non positivity of the scalar 
product raises substantial questions on selfadjointness and existence and uniqueness of time 
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evolution. In HS09] such problems have been treated on a slightly different version of 
the model, providing a general framework for existence and uniqueness of the Heisenbcrg 
time-evolution within a Hilbert-space formulation. 

The same methods could be used to discuss the evolution operators as unbounded op- 
erators on the Hilbert space, preserving the indefinite scalar product on a dense domain. 
Since, however, the Hamiltonian (|29[) is quadratic in the e.m. field variables and only in- 
volves a commutative algebra in the charged particle variables, a more pragmatic approach 
is enough for our purposes, with time-evolution operators defined in terms of Weyl expo- 
nentials of fields, introduced starting from their algebraic relations, on a suitable invariant 
vector space. 

Uniqueness of the solution will be obtained by observing that our space can be identified 
with a dense domain of the Hilbert space introduced in HS09J, that our time-evolution 
operators, although unbounded, are continuous and differentiable on Y in the Hilbert 
strong topology and that any group of isometries of "f , differentiable in the strong sense 
and with derivative given by (|29[) . coincides with them. 

The space Sf is defined as follows. Let srf ex™' be the *-algebra generated by the pho- 
ton canonical variables and by variables (Weyl operators in momentum space) W ( g , h ) , 
indexed by four- vector real- valued functions in L 2 ( K 3 ) , fulfilling 

W(g, h)* = W(-g,-h) , W(g , h)* W(g , h) = 1, (30) 

W(g, h) W(l,m) = exp(i((g,m)-(h,l))) W(l,m) W(g, h) , (31) 

[a(J),W(g, h)] =^= (f,n) W (g , h) , n = g + ih. (32) 
V 2 

In the above formulae the symbol * stands for the algebra involution and 

(f,g) = (f°,g°) - U\9 1 )- (33) 

i 

The Fock functional is characterized by the following expectations over stf ex™' 

u F (a(fi) a*(/ 2 )) =-</i,/ 2 >. (34) 

w F {W(g,h))= exp(±((g,g) + (h,h))). (35) 

In fact, expectations over monomials of a and a* can be expressed in terms of (|34[) with 
the aid of Wick's theorem |Wic50j , while those over monomials of W are identified by (1331) 
up to a phase factor, given by (13"T1) . and the other expectations follow from (1321 . 

The Gelfand-Naimark-Segal (GNS) construction |Nai59] on a *-algebra A and a linear 
(non-positive) functional ui proceeds as in the positive case, resulting in a non-degenerate 
indefinite vector space, a representation of A on it, with the * operation represented by 
the indefinite-space adjoint t , and expectations over a cyclic vector representing u> . 

The space Sf is obtained via the GNS construction over u> F ( srf ^ t m ' ) and its indefinite 
inner product is denoted by (■,■)■ In the sequel, ~f is regarded as a topological space with 
the weak topology r w , defined by the family of seminorms p (x) = \(y ,x)\ , with y e "V . 
The domain Y Q is given by 

r = J^(R 3 ) <g> Sf , (36) 

with Sf o obtained by performing the same construction described above for Sf , for canonical 
and Weyl operators with test functions in ( M 3 ) . Since no confusion should arise, the 
indefinite inner product of Sf is again denoted by (.,.). Sf is weakly dense in Sf and 
"V q is weakly dense in V by density of the Schwartz space in L 2 and Schwartz's inequality, 
applied to the explicit expression of the inner product. 

The representation of the Weyl operators defines the corresponding exponentials of the 
creation and annihilation operators on Sf ; such exponentials are also given by their series, 
which converge weakly on Sf , a fact which is not needed in the following. The above 
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construction is suggested by the fact that the formal time-evolution operator defined by the 
Hamiltonian (p?5)) contains exponentials of the canonical variables of the soft-photon field. 

We recall that indefinite-metric spaces obtained via a GNS procedure, starting from 
a non-positive functional and a *-algebra s$ ', are neither complete nor do they admit a 
unique completion. In some generality, complete spaces can be introduced as Hilbert-space 
completions of vector spaces obtained through GNS constructions. While in general the 
Hilbert structure can be very relevant for the existence and control of limits (the role of 
Hilbert space structures in the formulation of models of indefinite metric quantum field 
theories has been discussed at length in [MS 80] ). in the case of simple models, only involving 
polynomials and exponentials of fields, the use of a vector space looks simpler and even more 
intrinsic; a strong topology will only be needed for the formulation and control of uniqueness 
of time-evolution operators. 

The standard positive scalar product on Sf , obtained through the change of the sign of 
the summation over the space indices on the r.h.s. of (f55|) . gives rise to the Hilbert space 

introduced in }HS09j . 

The free and full dynamics of the model are determined by isometries of V leaving Y Q 
invariant and differentiable on "f in the strong topology of Jff , with derivative respec- 
tively given by H and H ^ PFBR ^ 

Notice that isometries U of a non-degenerate indefinite space 3f are determined by their 
restriction to a weakly dense subspace; in fact, by non degeneracy, U x , x £ 3? , is deter- 
mined by 

(y,Ux) = (If- 1 y, x) = lim ( t/ ~ 1 y, x n ) = lim ( y , U x n ) , 

n n 

with y G £?f, x n in a weakly dense subspace. 

The uniqueness of the evolution operators is a consequence of the following observa- 
tion, only requiring the existence of a positive scalar product (. , .) majorizing the indefinite 
product, \(y , x) \ < {x , x) 1 / 2 (y , y) 1/2 : 

Lemma 1. If two one-parameter families of isometries U(a) and V (a) of a vector space 
V q , endowed with a non- degenerate indefinite inner product ( . , . ) which is majorized by a 
positive scalar product, are differentiable on V in the corresponding strong Hilbert topology 
with the same derivative — i H (a) , then they coincide. 

Proof. Since strong differentiability implies strong continuity and hermiticity of H ( a ) on 
V , one has, V x , y G V , 

/- (x, V(-a) U(a) y) = -A (V(a) x, U (a) y) 
da da 

= -i (V(a) x, H(a) U (a) y) + i (H (a) V (a) x , U (a) y) = 0. 

□ 

The dynamics of the model in the interaction representation is determined by an operator 
U j : y — > Y such that the spectral component U I v , obtained from the decomposition 
of U j with respect to the joint spectrum of the components of v , is (.,.)- isometric for 
(almost) all v and leaves invariant a T w - dense subspace of "V , on which it is strongly 
differentiable with a time-derivative satisfying (I13|) . By employing formula (|14[) . one obtains 
a formal solution of the form 

C/,»(t) = exp(-ie (a' (/«(*)) +a(7<(*)))), (37) 

mKt) = I^£Lsl^±. (38) 

The operators (1571) fulfill (fT5j) in the strong topology, on f , which is left invariant. 
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Proceeding as for the PFB model, we introduce the adiabatic mean and the renormal- 
ization counterterm, 



H^ ] R = m + - m v 2 + iJ e ' m - + e v ■ A (p, x = 0) e" e 1*1 

3 9 _ 9 « I t I 1 



^ e 2 m z e^ 2 £ 1*1 ^ m + - m v 2 + em - + ffg « , (39) 

and determine the corresponding evolution operator U \ £ \ ( t ) in the interaction picture. 
For positive times, U \ e \ (t) is of the form 

U<r e < \(t) = hi'Ht) exp(- ?e (at(/( e )(0) + a(7- (e) (0))), (40) 

with 

/| i(k,t = P L— M ^— - — ■ (41) 

Considering that the above interaction Hamiltonian is Lorentz covariant only up to the 
second order in the four-velocity of the particle, it is not surprising that the above mass 
renormalization term, of the same form as that occurring in QED , only allows for the 
construction of the Moller operators up to higher (fourth order) terms in the four-velocity. 
We drop therefore the fourth order terms in h z and consider the (interaction representation) 
evolution defined by (|iP]). (|4Tj) . with 

h { z e) (t) = exp (- d (e) (*)) exp (-3ie 2 mz 6 * - 1 ) . (42) 

2 8 e 

For t < , one has to replace e by — e in the expressions above. 

The asymptotic limits of such evolution operators are given by the following 

Proposition 2. The large-time limits and the adiabatic limit of the evolution operator |^0| ) 
exist on "V and define Moller operators as isometries of "f : 



fl ( + X) = t,„ - lim lim U\ e \(-t) 

= h TiZ exp( ie [at (/^)+ a (/,-)]), (43) 
h T>2 =lim n lim /i^) (0 = 1 + O (v 2 ) , (44) 

/nk) = 4^-L. (45) 

Proof. By pT1) - ([55|) . the convergence of the coherence functions (|41|1 in L 2 implies the weak 
convergence of the corresponding expectations on Sf of polynomials and exponentials of 
the smeared photon variables. f2 £ A ' are invertible and preserve the inner product ( . , . ) 
as a consequence of the GNS representation of (131)1) . ([3~Tj) . They define therefore isometries 
of Y , uniquely determined by their restrictions to Y . □ 

As before, one can define a scattering operator 

S (PFB R) = u( e) {t) w { t>) (46) 

taking W as an isometric operator on L 2 ( R 3 ) (g> £f , with the same role as in (|2"51) ; in 
the comparison with the diagrammatic expansion, its matrix elements are interpreted as the 
non-infrared contributions to the corresponding scattering process. 

The expansion in powers of the electric charge of suitable matrix elements of the Moller 
operators of the PFBR model reproduces qualitatively the infrared contributions of Dyson's 
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power series, with modifications only due to the dipole approximation and to the non- 
relativistic limit. For instance, consider the transition amplitude for the scattering between 
two single-particle states ip v and tj) v , , without external (massive) photons, 

= <o ( _ A) v ,* F , n { + x \* F ) r v , tV , (47) 

where fl ^_ A \ is obtained from the decomposition of ^_ A with respect to the joint spectrum 
of the components of v . The radiative soft-photon corrections to the basic process are 
reproduced by the term in brackets in the second line of (|4"T|) . The exponentiation of the 
low-energy radiative corrections is thus a consequence of (|14l) and the compatibility of the 
dipole approximation with the renormalization procedure, with the limitation to the second 
order in the particle four-velocity discussed above, is also displayed. 

However, the transition amplitude (|47[) is not equal to the corresponding expression in 
the Coulomb gauge, given by ( S ^ PFB ^ ) v , v , namely it is not gauge invariant. For clarity 
of exposition, we point out that statements about the property of gauge-independence are 
not referred in this work to the concept of gauge invariance in its broad sense. Rather, since 
we employ the interaction representation and the adiabatic approximation, the asymptotic 
states obey a free evolution; the notion of gauge-independence of the transition amplitudes 
between physical scattering states in our models thus coincides with that of a free theory, 
namely with "the cancellation of the contributions from unphysical polarizations" to such 
amplitudes. 

A first indication that problems may arise comes from the fact that the dipole ap- 
proximation prevents the local conservation of the electric charge in a gauge involving four 
independent photon degrees of freedom; as a consequence, the standard argument explaining 
the cancellation of the contributions from longitudinal and scalar photons to the Feynman 
amplitudes no longer applies, since it requires the free-field character of d ■ A , which in 
turn relies on the continuity equation. The occurrence of infrared effects, due to the lack of 
(local) charge-conservation, can be explicitly checked; a simple calculation yields in fact 

(n^,*,,^*,) = exp(-^ |v'- v| 2 / ^* p 2 (k)) 

+ exp(-^ |v'- v| 2 / ^ p 2 (k)). (48) 

An extra factor 3/2 thus arises in the first exponent on the r.h.s. of (14"81 with respect 
to the exponent on the second line, which is obtained by evaluating the scalar product 
(fl_ tr V ' *f ' ^ + A .\r, v ) involving the Moller operators of the PFB model; con- 
sequently one has a different behaviour for the vanishing of (|47l) in the limit A — > , with 
respect to the corresponding Coulomb-gauge amplitude. The extra factor arises basically 
from the fact that the space components of the coherence functions (|45l) are of first order in 
the velocity, while only two components appear in (|27|). 

The relevant point is that the lack of gauge-invariance of (|4T|) depends neither on the 
dropping of the second order term in the interaction Hamiltonian nor on the approximation 
adopted in the computation of the evolution operator (|40[) . but is rather a difficulty common 
to all four-vector gauge quantizations in the presence of a dipole approximation. 

In order to trace back the origin of such effects, it suffices to notice that the expansion of 
Dyson's S- matrix elements for, say, the scattering of an electron by a potential contains, in 
the non-relativistic limit, powers of the scalar photon term (p ' / (p ' ■ k ) — p / (p ■ k ) ) 2 ~ 
oj ^ 2 [ ( v ' — v) ■ k] 2 , which do not appear in the presence of the dipole approximation; 
the residual contributions to the corresponding scattering amplitudes are thus due to the 
longitudinal photons. 

In HS09 , it has been pointed out that the dipole approximation prevents a straightfor- 
ward application of the Gupta-Bleuler procedure, since the d ■ A field no longer obeys a free 
evolution. As stressed in the above discussion, the difficulties with such an approximation 
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are indeed not limited to the implementation of the Gupta-Bleuler condition, but are even 
more basic; spurious gauge- dependent contributions, also relevant in the infrared regime, 
are in fact introduced into the Moller operators, thus not allowing to reproduce substantial 
features of the perturbativc treatment of soft-photon effects in a local and covariant gauge. 
The discrepancy with the perturbative results extends |PhDth| to the explicit expressions 
of the infrared- finite inclusive cross-sections |Wei95[ IJR76] . 

As it will be shown in the forthcoming Section, such difficulties do not arise in models 
based on a Bloch-Nordsieck approximation. 

2 Bloch-Nordsieck Models 

In the present Section we introduce hamiltonian models based on an approximation first 
devised by Bloch and Nordsieck [BN37 , which turns out to amount to a first-order expan- 
sion around a fixed four-momentum of each charged particle, with respect to the energy- 
momentum transfer. Since this fact is not transparent in the original treatment, a brief 
discussion may be useful. 

Consider the one-particle Dirac Hamiltonian with minimal coupling, 

H = a-(p-eA) + (3m+eA° = H D - e a ■ A + e A° , (49) 

and an eigenstate ip + p ( x ) = e~ l p ' x u r ( p ) of H D with momentum p and positive 
energy E , u r ( p ) being a spinor with helicity r . 

Let u r (p)=M r (p ) + O ( p — p o ) ; by the algebraic relations for Dirac's matrices one 
finds 

ff fl ^p(i)=[v'P + Vl- v 2 m] i> 0>p (x) + O(p-p ), (50) 
$o,p( x ) = e- ip ' x u r (p ) , v = £ Po /p . 

The v - dependent terms on the r.h.s of (|50p could also be obtained by formally replacing the 
matrices a and (3 in H D respectively by the (diagonal in the spinor indices) matrices v and 
V 1 — v 2 . Although this result may seem to rely on the linearity of H D with respect to the 
a matrices, it is indeed more general; for instance, it would also be obtained by performing 
a similar expansion on the eigenvalue equation for the Klein-Gordon Hamiltonian. 

According to the above discussion, we introduce the models defined by the Hamiltonians, 
respectively in the Coulomb gauge and in the FGB gauge, 

H f = p • v + H*- m tr - e v • A tr (p, x) = H% + H g t , (51) 

= p • v + Hf m - + e t> ■ j4 (p , x) = H q + H f nt , (52) 

with v = ( 1 , v ) , v a triple of self-adjoint operators in a Hilbert space, to be identified 
as the observable corresponding to the asymptotic velocity of the particle. They commute 
with the Weyl algebra s4 h generated by the canonical variables of the electron and with 
the polynomial algebras generated in the Feynman gauge by the photon canonical variables. 
Due to the appearance of the operators x , the above Hamiltonians do commute with space 
translations. The e.m. potentials occurring in (|51[) . (|52l) will be interpreted as describing 
soft-photon degrees of freedom. 

First we consider the model formulated in the Coulomb gauge. With the same notation 
as in Section [TJ the Hilbert space is 

,¥f = L 2 (M 3 ) <gs L 2 (R 3 d 3 v) <g) J? = Jd 3 vJf v , 

v is a multiple of the identity on Jf? v and the model can be studied at fixed v , with 
Hamiltonian H , v ' given by (ISTl) . 

Proceeding in a similar way as in the proof of the self-adjointness of the Pauli-Fierz 
Hamiltonian, one finds that H \^ is e.s.a. on D C J^ v ; existence and uniqueness of 
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the dynamics thus follow by Stone's theorem. The equations of the motion in the interaction 
representation, with the aid of (fTi|) . give 

&i,tr,x(t) = c x (t) exp (ie (a t * r (/ vx (0) + 0*r(7vx (*)))) , (53) 



c i ( t) = exp ( d 1 ( t ) ) , 

, / n /" g? 3 £; _ , . , . . sin v ■ k t . 

d x (t = / p 2 (k) (t — , 

J u) k v • k v ■ k 

n ( lr\ pi v ■ k t _ i 

'»■ ■'"■"^'■'•w— ' (M > 

With the same motivations and following the same treatment as in the first Section, we 
introduce the adiabatic renormalized Hamiltonian 

= H { Q v) +HiV e-^' + e 2 z x (v) v 2 e' 2 ^l 

= + *]>M«) , (55) 

1 f d 3 k p 2 (k) 
3 J uj k v ■ k 

Kato's result on time-dependent Hamiltonians applies as in Section [T] and, using (|14l) . one 
obtains the corresponding evolution operator in the interaction representation; for positive 
times 

*j, e l,A(*) = h ilHt) ex P ( ie (a* tr {f&(t)) + a^jjx (*))))> (57) 
ft(0(t) = exp(^^ rfi £) (0) exp( ?e 2 z 1 ( W ) v 2 ) , (58) 

di°(*) = - / U ^2 + 2, /5 2 (k)(e"- sin.-fc* 

J ^ k ((« • fc) 2 + € 2 ) 

+ ^ ie- 2 ** -1)), (59) 

n ( p i v ■ k t — e t _ 1 

/i e s ) x (k,o= e - lk x -£p= v 'M k ) — 1 i - (6°) 

v /2a; k i v ■ k — e 

The results for t < are obtained by replacing e with — e in the expressions above. One 
can then prove the following 

Proposition 3. The evolution operator in the interaction representation corresponding to 
the Hamiltonian \55}) , with the counterterm given by i56}) , admits asymptotic limits, yielding 
the Moller operators in the Coulomb gauge 



( A ) — o _ Hrv, Hrv, CO/ ( £ ) 



± , tr 



lim lim &7 e > x (-t) 
= exp(-ie^ [a* (/ v sx ) + a s (/ v , x )]) , (61) 



k-x P ( k ) * v • e s (k; 



/ vsx (k) = x -^=L ■ (62) 

^/2w k v ■ k 



The t — > ± oo and e — > limits in S61\) exists in the strong topology. 

Proof. The control of the limits is very similar to the case discussed in Proposition [1] □ 



13 



Scattering operators can be introduced as before in terms of operators acting on the "asymp- 
totic velocity Hilbert space" L 2 ( K 3 d 3 v ) . 

We shall now turn to the analysis of the model defined by the Hamiltonian ([52)l . which 
will be also referred to as four- vector BN model. Concerning the problems posed by the 
absence of a positive scalar product, we adopt the same choices as in Section [TJ following the 
same steps with identical results. The construction of the space of the model demands some 
modifications due to the effects of translation invariance, which requires the introduction 
of e.m. canonical and Weyl operators, fulfilling equations of the same form as (|30p - (|3"!2")) . 
smeared with test functions of the form 

<7 x (k) = cos(k-x) /(k), /i x (k) = - sin(k-x) /(k) , (63) 

with / ^ real valued and square integrable. 

The GNS construction over such an algebra and the Fock vacuum yields an indefinite 
space Sf BN and a weakly dense domain & B n ; the tensor product with the same particle 
space as for the Coulomb-gauge version of the model gives the indefinite-metric space 

Y BN = L 2 (R 3 ) ® L 2 {R 3 d 3 v) ® <S BN 

and the weakly dense subspace 

y BN = 2> (R 3 ) ® y (r 3 ) ® sr 

y BN and i /BN will be decomposed on the spectrum of v as above. 

Isometric evolution operators on "f BN are constructed as in Section [T] and are unique 
in the same sense. As for the Coulomb gauge, the model can be studied at fixed v . The 
evolution operator in the interaction representation is given by 

*M*)= c 3 (t) exp(-ze (at(/ ux (0) + a(T vx (t)))), (64) 

i c v 

c 2 ( t ) = exp ( d 1 {t)), 

J 2 lu k i v ■ k 

In order to construct Moller operators, we introduce the renormalized Hamiltonian and the 
adiabatic factor, 

Hi v) R ie) = p v +H Q e -J +e v ■ A(p,x) e - £ "l 

-e 2 z 2 (v) v 2 e- 2e l*l = ffSrU^^, (65) 



1 [ d 3 k p 2 (k) 3 

^ 2 ( « ) = o / r~ = o z i ( u ) ■ 66 ) 

2 7 w k V ■ K 2 

For t > , the corresponding evolution operator in the interaction representation is 

W^\{t) = hl'Ht) exp(ie(at(/ w (0(t)) + o(7ix (*)))). (67) 



e -2 6 t _ i 

x exp ( - i e 2 z 2 {v) v 2 — ) , (68) 



i v ■ k t — e t 



The results holding for t < are obtained by replacing e by — e in the expressions above. 
The asymptotic limits are controlled by the following 
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Proposition 4. The large-time limits and the adiabatic limit of the evolution operator |6'7| ), 
defining the Moller operators of the four-vector BN model, exist in the weak topology and 
define isometries in Y BN , given by 



( ± ] = t w - lim lim (-t) = t w - lim n { ± X) e 

= exp(*e [at(/_ x ) +a(J v _)]), (70) 



Proof. The proof is the same as that of Proposition [21 □ 

Scattering operators for the four-vector BN model, interpreted as the result of soft photon 
corrections to an isometric scattering operator W in the velocity space, can be introduced 
by 

SI FGB) = r w - lim lim ^}%(t) W V}' t {(t') 

£ — V t , t ' — ¥ +00 ' 

s ■ (72) 

More generally, for the comparison of (|72[) with the Feynman-Dyson expansion, the isometry 
W can be taken as acting on the entire space "V BN and interpreted as yielding an infrared- 
finite S - matrix, obtained by factoring out from the scattering amplitudes the contributions 
given by the Moller operators of the model. 

In the remainder of this Section, we prove that the infrared diagrammatic of QED 
is reproduced with the help of the operators (|70p . Regarding the infrared approximations 
and results within the perturbative-theoretic framework, we shall refer in the sequel to the 
streamlined treatment of |Wei951 IJR76] , while a more detailed analysis can be found in the 
classic work of Yennie, Frautschi and Suura [YFS61] . 

In what follows, we introduce an energy scale, say A , conventionally dividing the soft and 
the hard photons; for photon canonical operator-valued distributions smeared by solutions 
of the wave equation (with mass A ) with energy below and above A , we shall denote the 
corresponding subspaces of & respectively by & so t t and & hard . 

For defmiteness, we consider as a basic process a — > f3 the scattering of an electron 
by a potential and suppose that the incoming (outgoing) particle is described by a state of 
definite momentum, with four- velocity v ( v ' ) . The process is supposed not to involve low- 
energy photons, while it may involve hard photons. Let rj ™ m (?7°"m ) be the state of the 
incoming (outgoing) e.m. field; under the above assumptions, rjl n m° ut ^ = ^ F ® 7 l n m° ut ^ , 
with * F the vacuum vector of # ' soft and 7™ m .(7e.m.) belonging to & hard - 

The isometry W in (|72|) is interpreted as the hard-photon part of the scattering operator, 
hence it is supposed to act as the identity operator on & so t t ; moreover, according to the 
above interpretation of the models, the same property is assumed to hold for the restrictions 
of j_ A ' to & hard . The transition amplitude for the basic process a — > /3 is therefore 

Tr pa = {i> v ,®l° Ut m.^{i> V ®lTm.)), (73) 

with if) v improper vectors describing a charged particle of four- velocity v , associated to 
the spectral resolution of the operators v . 

The matrix element for the same process, including the contributions due to the Moller 
operators, is 

(si FGB) ) Pa = <n ( _ A) (w® vTD, * « ( + ) (^,®'?r m .)>- (74) 

The first result is given by: 

Proposition 5. The soft-photon radiative corrections to the process a — ¥ (3 are reproduced 
by the following contributions to the transition amplitude {7$ : 

(fiL*'* <* F , * F ) exp(e 2 [a(/„, x ),at (/_)]). (75) 
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Proof. The proof is based on direct calculations. 

The last term of (1751) is related to "virtual soft photons" emitted from either the incoming 
or the outgoing external fermion leg and absorbed from the other, and is equal to 

The first (second) expression in brackets in ([75)1 is related to the electron wave-function 
renormalization relative to the outgoing (incoming) electron line; it is given by 

X 2 



--(-T 3 (.V) '■-'»■ (77) 

with Z 2 x ( v ' ) the infrared part of the wave- function renormalization constant relative 

to the outgoing leg and £ ^ \ ( u ' ) the self-energy insertion on the same line, the electron 
mass being renormalized: 

"'»*■'>(''>- e2 °''/^- (1 -t ( vV +0(g2) " 

with k = k / a; k . Such a term is a residual infrared contribution, arising whenever "on- 
shell" mass-renormalization conditions are imposed |JR76] . 

Summing (176[) and (1771) . the standard result for the soft-photon radiative corrections 
((13. 2.4), (13.2.5) in |Wei95j i is recovered. □ 

We can now examine the effects due to the emission of low-energy radiation. For definite- 
ness, we consider the contributions due to the emission of n soft photons from the outgoing 
fermion line. Let jr 1 "t " denote the vector describing the n - photon state with polar- 
ization indices fi fi n and four- momenta k k n . It is straightforward to establish the 
following: 



Proposition 6. The overall transition amplitude for the process a — > (3, with emission of 

' k 1 "' k 



n soft photons described by the state ^ f 1 "u n , is reproduced by the matrix element 



( v> v ■ ® * t\ :X B ® t °. ™- ' 5 He) ( ^ t, ® * ® 7 r TO . ) ) ■ (78) 

Proof. A simple calculation yields 



2 w k - « ' • & j - * e 



(79) 



Recalling ([TO]), ([72} and (J7JJI and employing (JTHJ) one gets 

( *v ® ® iri . (>„ ® *f ® 7i n m .)> 



, q (FOB) x n n V^jJ ST 

v A , (e ) 11 3 = 1 / = 



v 1 = v,v 2 = v',r] 2 =l=-r] 1 , 
in agreement with the perturbative expression ((13.3.1) in [Wei95]). □ 
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The infrared phases occurring in the transition amplitude for a process described by the 
sum of one-particle Hamiltonians, with at least two charged particles in cither the initial or 
the final state, can also be recovered by explicit calculations. For brevity we do not report 
the details. 

Finally we remark that the Bloch-Nordsieck expansion preserves the free-held character 
of d ■ A and consequently the gauge-invariance of the transition amplitudes, in the sense 
that terms from unphysical photon polarizations cancel. 

Outlook 

By exploiting the control of the soft photon contributions to the Feynman-Dyson expan- 
sion of QED , gained through the four-vector Bloch-Nordsieck model, the recipes leading to 
infrarcd-hnitc inclusive cross-sections can also be formulated and discussed. In particular, 
it is possible to enlighten the non-perturbative limitations of such recipes, due to the su- 
perselection of particle momenta, taking advantage of the fact that within the hamiltonian 
approach here developed issues connected with an order-by-order diagrammatic treatment 
are avoided. We plan to report on these problems in a future paper. 
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